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New moderate-resolution mid-infrared spectroscopy at 10 /im of 27 infrared galax- 
ies is presented. The galaxies have been chosen from three 60 /xm selected and one 12 
/itm selected complete flux-limited catalogs of galaxies; 17 of these sources have Ltr(8- 
1000/xm) > 5 x 10 11 L Q . A high-resolution spectrum of the source Arp 299B1 is also 
presented. Combining these new results with previously published results, a nearly 
complete 60 fim selected flux limited subsample, with LiR(8-1000rmi) > 1.6 x 10 11 
L Q of 25 galaxies is defined. Within this subsample, it is found that the domi- 
nant power source of infrared galaxies in the luminosity range 1.6 x 10 11 < Lm,(8- 
1000/xm) < 10 12 L/0 is massive star formation based on the clear presence of the 
11.3 fim aromatic hydrocarbon emission feature in the majority of the spectra. Three 
of the five ultraluminous infrared galaxies (LiR,(8-1000/im) > 10 12 L Q ) within this 
subsample show evidence that an active galactic nucleus provides an energetically im- 
portant power source based on the detection of silicate absorption in their mid-infrared 
spectra. The physical basis of a possible anti-correlation between the 11.3 /xm feature 
equivalent width and infrared light to molecular gas mass ratios is discussed. 
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1 INTRODUCTION 



This work is aimed at the question of what provides the 
dominant power source of ultraluminous infrared galaxies, 
that is galaxies with luminosities in the wavelength range 
8-1000 fim in excess of 10 12 L (H =75 km s _1 Mpc" 1 ). 
This question is of interest because it is possible that the 
power source of these galaxies is qualitatively different from 
the dominant power source of most infrared galaxies of lower 
luminosity. The observational approach of this work is the 
use of moderate resolution 8-13 /im spectroscopy. 

In very broad strokes, the context of the debate over the 
question addressed here can be summarized as an attempt, 
on the one hand, to explain ultraluminous infrared galaxies 
as more powerful versions of starburst galaxies, that is as 
super starbursts (Joseph & Wright 1985), where massive 
star formation provides the power source, and on the other 
hand, to demonstrate that ultraluminous infrared galaxies 
are the precursors of quasars (Sanders et al. 1988) and are 
powered by active galactic nuclei (AGNs). 

The question has remained open to serious debate for 
essentially one reason: dust. Infrared galaxies by their nature 
radiate the bulk of their energy through the degradation of 
the optical and UV photons emitted by the primary power 
source into infrared photons emitted by dust heated to a 
temperature of ~50 K. This has the effect of erasing most 



of the commonly accepted features used to distinguish the 
characteristics of starbursts and AGNs. 

Optical and, to a lesser degree, near-infrared studies 
are particularly susceptible to this difficulty since spectro- 
scopic signatures of AGNs in ultraluminous infrared galax- 
ies should only betray sources that are relatively unobscured 
and therefore might not contribute substantially to the heat- 
ing of the dust that is the final source of most of the ob- 
served emission. However, Veilleux et al. (1995) have shown 
in a large sample of infrared galaxies that the fraction of 
infrared galaxies with AGN-like optical spectra increases as 
a function of increasing luminosity, providing important cir- 
cumstantial evidence that AGN-like activity is associated 
with ultraluminous infrared galaxies. The strength of this 
evidence might be considered to be weakened, however, by 
the finding of Goldader et al. (1995) that the ratio of Br7 
to infrared luminosity decreases as a function of luminos- 
ity, suggesting that dust obscuration plays an even greater 
role in ultraluminous infrared galaxies than in luminous in- 
frared galaxies. This suggests the perverse possibility that 
energetically dominant starbursts in ultraluminous infrared 
galaxies are so obscured in the optical that only incidental 
unobscured AGNs are detected. 

Radio observations are unaffected by dust obscuration 
and can provide high spatial resolution. However, the pri- 
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marily synchrotron radiation that is detected is only indi- 
rectly linked to the power source of infrared galaxies. In star- 
bursts, it is thought that the radio emission arises from elec- 
trons accelerated in supernovae remnants, while in AGN, rel- 
ativistic jets are strong radio emitters. Condon et al. (1991; 
hereafter CHYT) have employed the high spatial resolution 
provided by the Very Large Array (VLA) to show that the 
radio sizes of a large fraction of ultraluminous infrared galax- 
ies are comparable to the minimum size required for a black- 
body that reproduces the 25, 60, and 100 fim flux densities 
of these sources. From this they infer that starbursts rather 
than AGNs may power these sources. At even higher spa- 
tial resolution provided by very long baseline interferometry 
(VLBI), Lonsdale, Smith & Lonsdale (1993) find evidence 
for radio structure on the scale of milliarcseconds, some com- 
ponents of which have brightness temperatures well in ex- 
cess of what could be expected in starburst-related emission. 
They deduce from this that AGNs are present in five of-eight 
of the ultraluminous infrared galaxies in their sampleEI. 

In this work, the dust emission itself is examined for 
spectral clues that can reveal the nature of the power source 
in infrared galaxies. In a sample of 60 galaxies, Roche et 
al. (1991; hereafter RASW) have shown that their 8-13 
/jm spectra can be classified into three types: those with 
prominent polycyclic aromatic hydrocarbon (PAH) features 
generally show Hll region-like optical spectra, those with 
flat featureless spectra have Seyfert 1 or quasar-like opti- 
cal spectra, and those rare galaxies with silicate absorption 
features usually have Seyfert 2 optical spectra. Dudley & 
Wynn- Williams (1997) have extended the work on galax- 
ies with very deep silicate features to show that the size 
of the power source can be constrained to be smaller than 
a few parsecs, implying that such sources are powered by 
deeply embedded AGNs. Here, the 8-13 /xm spectra of a 
well-defined sample of 25 luminous and ultraluminous in- 
frared galaxies are examined in this context. 

It should be noted that all of the sources for which new 
spectral observations are reported here can be found in one 
or more of four IRAS-based catalogs: the Bright Galaxy Sur- 
vey (Soifer et al. 1989), the Extended Bright Galaxy Survey 
(Sanders et al. 1995), the Two Jansky Survey (Strauss et 
al. 1992), or the Extended 12 pm Survey (Rush, Makan & 
Spinogio 1993). 

Observations and data reduction procedures are de- 
scribed in Section 2. Section 3 presents the new spectral 
observations. The classification of the spectra is discussed 
in Section 4.1. The justification for the classification of the 
ultraluminous galaxies that fall within the later subsample 
is discussed in some detail in Section 4.2. In Section 4.3 
a subsample of galaxies is defined and an apparent trend 
with luminosity within it is examined. The possibility of an 
anti-correlation between the equivalent width of the 11.3 
/im PAH feature and the ratio of infrared light to molecu- 
lar gas mass is discussed in Section 4.4, and conclusions are 
drawn in Section 5. 



* Smith, Lonsdale & Lonsdale (1998) suggest that radio super- 
nova remnants may be responsible for some of the structure they 
detect 



2 OBSERVATIONS AND DATA REDUCTION 

The new observations reported here were obtained over the 
course of 5 years, from 1991 through 1995, using the 10 and 
20 /im Cooled Grating Spectrograph (CGS3) mounted at 
the Cassegrain focus of the United Kingdom Infrared Tele- 
scope (UKIRT) on Mauna Kea in Hawaii. Observations in 
1991 were obtained by G. Wynn- Williams and J. Goldader, 
in 1992 by G. Wynn- Williams and C. Dudley, and in 1993- 
1995 by C. Dudley. For parts of each observing run a UKIRT 
scientist (J. Davies, T. Geballe, or G. Wright) was present 
at the summit and participated in the observations. CGS3 
has been described by Cohen & Davies (1995). All but one 
of the spectra reported here were obtained using the low- 
resolution spectroscopy mode with A/ A A ~ 60 in a series of 
beam-switched, chopped observations with a chop frequency 
of ~10 Hz and a beam switch frequency of ~0.5 Hz. To fully 
sample each spectrum, two interlacing grating positions were 
observed with changes between the two positions occurring 
after 10-16 beam-switched pairs. Table 1 gives a log of ob- 
servations for the galaxies. In Table 1, column 1 is the object 
name, columns 2 and 3 give the aperture centre in RA and 
DEC, column 4 gives the position reference, columns 5 and 
6 give the object redshift and reference, column 7 gives the 
date of observation, column 8 gives the circular aperture size 
used in arcseconds (full width at 10 per cent power), column 
9 gives the chop throw in arcseconds, and column 10 gives 
the chop direction. 

Observations were calibrated in wavelength by means 
of high order spectral measurements of a Kr lamp in the 
mariner discussed by Hanner, Brooke & Tokunaga (1995). 
A master calibration spectrum was first selected and cali- 
brated in wavelength by flattening a 2x over-sampled Kr 
lamp spectrum (4 grating positions) with the spectrum of a 
soldering iron replacing the Kr lamp in the beam. Kr lamp 
spectra from each night's observations were then cross cor- 
related against the unflattened master spectrum to give a 
wavelength calibration for each night. Typical shifts from 
the nominal CGS3 wavelength calibration were less than a 
resolution element, and were never more than 2 resolution 
elements. 

Corrections for instrumental response, atmospheric 
transmission, and flux calibration were performed by divid- 
ing source spectra by the spectrum of a standard star with a 
known 10 fim flux density and an assumed blackbody tem- 
perature based on its spectral type, and then multiplying by 
the Planck curve of that temperature scaled to match the 
standard star in flux density. Observations of standard stars 
were performed with the same aperture, chop throw, posi- 
tion angle and frequency as the galaxy observations they cal- 
ibrate. The standard star observations were made to match 
in airmass so that no additional airmass corrections would 
be required. For the standard stars, the number of beam- 
switched pairs between grating shifts was typically 3-5. The 
stars have been chosen, where possible, to have a spectral 
type of K0 or earlier to avoid the effects of SiO absorption, 
H~ opacity, and circumstellar dust emission that can effect 
the spectral shape of late-type stars. Table 2 gives the HR 
number, 10.1 /im flux density, photometric reference, and 
assumed blackbody temperature of each standard star, in 
columns 4, 5, 6, and 7, respectively, for each of the observa- 
tions given in Table 1. In column 2 of Table 2, alternative 
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names are given for galaxies that are found in multiple cata- 
logs, while column 3 gives the date of each observation as in 
Table 1. It should be noted that an error of ~1000 K in the 
assumed blackbody temperature of the standard stars leads 
to a ~2 per cent error in flux calibration at either end of 
the passband 8-13 fim which is negligible for the purposes 
of this work, while for those stars of late spectral type, de- 
viations from the assumed blackbody can be ~10 per cent 
in the SiO band (Cohen et al. 1995). Systematic uncertain- 
ties in overall flux calibration of 15-20 per cent can arise 
due to telescope pointing errors with the stronger effects for 
the smaller apertures. For wavelengths shortward of 8.2 /jm, 
strong systematic effects can occur as a result of a variable 
water vapor column in the Earth's atmosphere. 

All spectra were regridded to a common wavelength 
scale using a Gaussian filter of a = 0.064 /im , which lead 
to a slight correlation between neighboring data points. 
Those spectra of the same source obtained on two or more 
nights were combined, weighted by integration time. Data 
not falling in the ranges 7.7-9.1 and 10.0-13.4 /im have been 
excluded due to poor atmospheric transmission. 

Spectra of Wynn- Williams et al.'s (1991) Arp 299B1 
were obtained on 1994 February 6 and 7 using CGS3 in 
high-resolution mode (A/AA ~ 250) using a 3'.'26 full width 
at 10 per cent power circular aperture. The spectra obtained 
on each night were the result of two sets of observations each 
consisting of two interlacing grating positions but with dif- 
ferent wavelength ranges to cover the spectral range 10-13.2 
/jm. The chop throw was 20", and the chop and beam switch 
PA were 45°. About 30 per cent of the first night's observa- 
tions were affected by clouds; these data were rejected by in- 
specting the change in the dispersion of the total signal. The 
Kr lamp was observed on each night to provide wavelength 
calibration. The strongest lines from each night agreed in 
wavelength to 3 per cent of a resolution element so that the 
Kr spectra of the two nights were combined, and a differ- 
ence in wavelength of 30 per cent of a resolution element 
was measured between the expected (based on the nominal 
calibration of CGS3) and the measured wavelengths of 5 Kr 
lines. An additive correction of 0.02 fim was applied to the 
observations, this being the interpolation of the 2.117 and 
2.191 /jm Kr lines measured at 6th order at the expected 
position of the galaxy [Ne II] line. Flux calibration was per- 
formed using HR 4067 with a 10.1 (jm flux density of 103 Jy 
and an assumed temperature of 4000 K on the first night, 
and HR 5340 using a 10.1 fim flux density of 738 Jy and an 
assumed temperature of 4400 K on the second night. No sig- 
nificant difference was found between the two nights' obser- 
vations, so they were combined according to their statistical 
weights. The difference in airmass between the source and 
standard observations was ~0.5 but monitoring of BS 4518 
on the first night at airmasses similar to the source obser- 
vations showed that the level of spectroscopic variation over 
the observed wavelength range as a function of airmass was 
negligible compared to the final signal-to-noise ratio. The 
overall photometric calibration is estimated to be good to 
about 25 per cent while the typical signal-to-noise ratio of 
an independent data point is about 10. 



3 RESULTS 

The intensity and equivalent width of the 11.3 fim PAH fea- 
ture have been measured directly for all of the sources by 
estimating the continuum in the rest wavelength (see Table 
1) ranges 10.5-10.9 and 11.7-12.1 fim and subtracting the 
linear interpolation of these two estimates from the spectra 
to measure the residual flux in the range 10.9-11.7 /im. Un- 
certainties in the feature intensities represent the combined 
uncertainties in the residual flux in this range but not the 
uncertainties in continuum subtraction. For the equivalent 
width estimates, uncertainties in the continuum level have 
been included. In most cases the continuum uncertainties 
are dominant in the equivalent width estimates, leading to 
indeterminate results in the case of limits, but in a few cases 
the uncertainty in the residual is the dominant component, 
and upper limits to these equivalent widths can be set. 

Since model fitting can lead to substantially cleaner 
estimates of feature fluxes (of order the signal-to-noise ra- 
tio in the peak spectral element being fit) estimates based 
on the model that is applied to these data to test for the 
presence of a silicate absorption feature are also given. The 
model consists of a power-law continuum and an additive 
PAH spectrum represented by the spectrum of the Orion 
Bar taken from Roche (1989). The fitting of this model 
was accomplished using a version of Bevington's (1969) pro- 
gram CURVEFIT supplied with IDL ver. 4.0 using three 
free parameters, the power-law continuum scaling and index 
(/„ = /3 x f~ a ), and the scaling of the Orion Bar spectrum. 
The data were fit in the rest wavelength ranges 10.0-10.9, 
11.0-11.4, and 11.7-12.6 /im using initial guesses based on 
linear fits to the logarithm of data binned in the ranges 10.0- 
10.5, 10.5-10.9, and 11.7-12.6 fim for the power-law scaling 
and index, and based on the continuum subtracted estimates 
of the flux in the 11.3 /im feature. The weighting of the 
data was statistical. As discussed recently by Goldader et 
al. (1995), the use of a power-law as a fitting function tends 
to throw the bulk of the errors in the fit parameters into 
the power-law scaling, so that the formal errors of the re- 
maining parameters are underestimated. Further, the model 
assumptions introduce further systematic uncertainties. In 
particular, by fitting the data in the range 11.2-12.6 fim, it 
is assumed that the ratio of 11.3 fim feature emission with 
respect to the underlying plateau emission that extends to 
12.9 fim is the same as that found at the position of peak 
of the 11.3 fim emission feature in the Orion Bar. This ra- 
tio, however, is know to vary spatially within the Orion Bar 
(Roche et al. 1989), and also within at least one external 
galaxy (Dudley & Wynn- Williams 1999). The consequences 
of this uncertainty are largest in spectra where the PAH 
emission is most powerful relative to the continuum. Given 
that the galaxy redshifts make it impossible to estimate the 
continuum beyond 13.1 fim (rest frame) in ~50 per cent of 
the sources, it has not been possible to quantify the im- 
portance of this effect in these data, regardless of obvious 
signal-to-noise issues. Due to these multiple and tangled un- 
certainties, both results of direct measurements and model 
fits with formal errors are presented, but it must be cau- 
tioned that the model parameter errors are seriously under- 
estimated. In particular, model 3 a upper limits to the 11.3 
fim intensity are roughly the same as measured 1 a upper 
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limits to the 11.3 /im intensity. The new spectra and model 
fits are presented in Fig. 1. 

Table 3 presents the model power-law index in column 
2, the measured (upper) and model (lower) 11.3 /im feature 
intensity in column 3, the 11.3 /jm feature equivalent width 
in column 4, the ratio of the flux density of the source in the 
rest wavelength range 8.0-8.4 /im to the model flux density 
in the same range where the errors represent the errors in 
the estimates of the source flux densities only. Measurements 
of fine structure lines are given in columns 6 and 7. For the 
measurements of fine structure lines the continuum was esti- 
mated by eye, and for detections of the 12.8 /im [Nell] line 
it must be cautioned that a significant amount of the in- 
tensity could be due to a blended PAH emission feature at 
12.7 (jm except in the case of Arp 299B1, where the spec- 
tral resolution is sufficient to distinguish the two features. 
Errors for detections and limits are 1 a and 3 a, respec- 
tively, except in the case of the model 11.3 fim equivalent 
width, which is reported in all cases to give a sense of the 
model continuum uncertainties. In the case of Arp 299B1, 
the fine structure [S iv] and [Ne II] results are derived from 
the high-resolution spectrum. The redshift of the [Nell] line 
(2900 km s -1 ) agrees within 1.4 a (1 a = 170 km s -1 based 
on Gaussian fitting) with the CO redshift measured by Sar- 
gent & Scoville (1991), while much more precise agreement 
between unpublished IRTF CSHELL Br7 data and the CO 
observations is known, so that some confidence in [Nell] 
identification is warranted. It should be noted that [Nell] 
is not detected in the low resolution 8-13 /im spectrum of 
Arp 299B1 but that the 3 a upper limit for this feature is 
2.1 xl0~ 15 W mT 2 , which is consistent with the intensity 
measured at high resolution. Further, the ratio of 1(11.3) to 
I([NeIl]) (1.7) is not unusual for galaxies for which the ratio 
has previously been measured (RASW). 

As is apparent in Table 3, while a wide range of power- 
law indices are necessary to describe the long-wavelength 
continuum, the ratio of 8.2 /im flux to model 8.2 /im flux is 
not significantly different from unity for a large fraction of 
the galaxies even in cases where the power-law index is large 
and a minimum near 9.7 /im is pronounced. As a result, it 
is not generally necessary to invoke silicate absorption to 
explain the observed spectra. 



4 DISCUSSION 

4.1 Spectral classification 

RASW have found that a three-type classification scheme 
is sufficient to classify the 8-13 /im spectra of most galax- 
ies. Their scheme consists of flat-featureless, silicate absorp- 
tion or emission, and PAH-type (their UIR) spectra. Fur- 
ther, there is a general correspondence between the 8-13 
/im spectrum of a galaxy and its optical spectral classifi- 
cation based on optical emission lines. The correspondence 
is that those galaxies with flat 8-13 /tm spectra are often 
broad-line AGNs, those with silicate absorption are fre- 
quently narrow-line AGNs, and those with PAH features 
typically have optical spectra with line ratios similar to 
Galactic H II regions. 

Parenthetically, the identification of the mid-infrared 
emission feature analysed here is not entirely secure. Other 



laboratory analogs, including quenched carbonaceous com- 
posites (Sakata & Wada 1989), hydrogenated amorphous 
carbon (Duley 1989), and anthracite (Ellis et al. 1994), have 
been proposed in addition to the PAH model (Puget & 
Leger 1989) that is adopted here for ease of expression. How- 
ever, the precise identification of the chemical source of these 
features is not particularly important for the present work. 
What is important is that the strength of one feature may 
be predicted to some extent from observations of another. 
In particular the strength of the 7.7 /im feature may be pre- 
dicted from observations of the 11.3 /im feature to within 
about a factor of 2 (Zavagno, Cox & Baluteau 1992). 

The model fitting performed at the long- wavelength end 
of the spectra can be used effectively to test for the pres- 
ence of strong silicate absorption. Two cases might occur. 
First, a PAH spectrum might undergo external extinction. 
In this case, two effects will contribute to the model under 
predicting the observed flux in the shorter wavelength half 
of a spectrum. The first effect is that the 11.3 /im feature 
will undergo more severe extinction than the region near 8 
/tm due to the spectral shape the silicate component of the 
interstellar extinction curve. This has the effect of reduc- 
ing the predicted PAH strength. The second effect is that 
the slope of the model continuum will have a larger spec- 
tral index than would otherwise be the case if no silicate 
absorption were present. The model would therefore under- 
predict the continuum underlying the 7.7 and 8.6 /im feature 
complex. Both of these effects will cause the model to under- 
predict the shorter wavelength observations. Second, a PAH 
spectrum may undergo no particular extinction, however, 
a second silicate-absorbed continuum component may also 
contribute to the spectrum. In this case, the model should 
predict the PAH component well. However, if the silicate ab- 
sorption is strong, the model spectral index will be increased 
due to the additional component, and the continuum in the 
short wavelength half will be underestimated so that the 
model will also under-predict the observed flux. 

The RASW classification scheme is applied to the 
present data in the following way: the ratio of the observed 
to model flux density between 8.0 and 8.4 is used as a guide 
to classification such that if this ratio is 3 a above 2.5, then 
the source is classified as silicate. This cut may somewhat 
underestimate the number of galaxies in which silicate ab- 
sorption is present compared to the approach taken in Dud- 
ley & Wynn- Williams (1999), where the cut is taken at 2.0. 
For the purposes of definite classification, under the work- 
ing hypothesis that all infrared galaxies are PAH-type unless 
convincingly shown to be otherwise, this more conservative 
value is appropriate. In cases where the 11.3 /im feature is 
not formally detected but only allowed (dot dashed curves 
in Fig. 1), and where a minimum is present near 9.7 /im, 
a question mark is appended to the PAH designation; how- 
ever, whether or not there is a formal detection of the 11.3 
/im feature in sources with obviously flat spectra, the clas- 
sification is given as flat. A mixed classification of PAH and 
silicate is given where PAH emission can explain ~30 per 
cent of the observed emission between 8.0 and 8.4 /im. Un- 
der this condition, the classification of Arp 220 by RASW 
as both PAH and silicate would be revised to just silicate. 
Arp 220 is discussed in detail in Section 4.2.5. 

Table 4 gives, for each source in column 1, the spectral 
classification based on the present data and on optical ob- 
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servations, along with 12-100 /im photometry, the infrared 
luminosity, and the star formation efficiency. In some cases 
the optical spectral classification has been inferred from the 
data of Veilleux et al. (1995), if some diagnostics were miss- 
ing. It is clear from Table 4 that the general correspondence 
between optical spectroscopy and 8-13 fim spectroscopy is 
quite similar to that found by RASW. Neglecting LINER- 
type galaxies, four exceptions where AGNs are identified op- 
tically but which are found to have PAH-type spectra in this 
work are Mkn 1034, IRAS 04154+1755, IRAS 05189-2524, 
and Zw 475.056. 

Mkn 1034 is a known Seyfert 1 galaxy; however its 8-13 
fim spectrum can be fit with a PAH plus continuum model. 
This may be consistent with the finding of Sopp & Alexan- 
der (1992) that only half of the 6-cm flux density measured 
in a 14'.' 5 beam is also contained within a l'/5 beam, sug- 
gesting that spatially extended star formation plays a role 
in this source. A similar explanation may be invoked for 
the presence of PAH emission in the spectra of the Seyfert 
2 galaxies IRAS 04154+1755 and Zw 475.056. Crawford et 
al. (1996) find IRAS 04154+1755 to be resolved at 6 cm 
using both 4-arcsec and 0.5-arcsec beams. In the map of 
Zw 475.056 presented by CHYT, low-level radio emission 
extends over 1.5-arcsec. IRAS 05189—2524 is discussed in 
detail in the next section; however it should be noted here 
that the presence of exceptions, in this sample, to the the 
correspondence between optical and 8-13 fim spectral clas- 
sification found by RASW should not be taken as evidence 
that the correspondence is merely coincidental. The sources 
presented by RASW were sometimes selected on the ba- 
sis of optical spectroscopy to examine the question of what 
kind of spectrum they would present at 8-13 fim. Hence, 
many of their sources were archetypical examples of their 
optical class. In the present sample, selection is based solely 
on infrared properties, so that mixed cases might well be 
expected. 

It is probably not coincidental that these four sources 
are among the sources in Table 3. with the smaller values 
of a suggesting that the AGN in these sources may con- 
tribute somewhat to the 8-13 fim continuum; strong 12 and 
25 /jm emission is known to select for AGNs that have un- 
ambiguous optical spectra (de Grijp et al. 1985; Spinoglio 
& Malkan, see also the AGN spectral energy distributions 
presented by Sanders et al. 1989). 

4.2 Comments on individual ultraluminous 
infrared galaxies 

The observed properties of the ultraluminous infrared galax- 
ies presented here show a surprising range. A closer look at 
these properties will be discussed in this section. 

4-2.1 A phenomenological description of IRAS 
05189-2524 

One of the greatest surprises in these data is the case of 
IRAS 05189—2524, which might have been expected to show 
a fiat or silicate absorbed 8-13 micron spectrum based on 
its optical AGN- like spectrum (Sanders et al. 1988). What 
is observed however is a spectrum that is well fit by a 
PAH plus power-law continuum model (see Fig. 1). Addi- 
tionally, the strength of the 11.3 fim emission relative to 



the 8-1000 fim emission is quite similar to known starburst 
galaxies so that it seems possible that star formation might 
be the primary power source of the galaxy. Here, a ten- 
tative decomposition of its spectral energy distribution is 
proposed. In Fig. 2, a comparison is made between IRAS 
05189—2524 and two known AGNs to suggest that while 
the 1-5 fim emission in IRAS 05189—2524 may be due to 
hot dust heated by an AGN, the 8-100 fim emission may be 
better explained as arising from star formation traced by the 
observed PAH emission. In Fig. 2 the power-law extrapola- 
tions of the dust continuum emission (solid lines) found in 
the 8-13 fim spectra are projected to 5 and 25 fim. For the 
two sources where the 8-13 micron emission is thought to be 
due to AGN-heated dust (Mkn 231 and IRAS 13349+2438) 
(Roche et al. 1983; Beichman et al. 1986), these extrapola- 
tions are reasonable predictors of the observed flux densities 
at both 5 and 25 ^m. For IRAS 05189-2524 the extrapola- 
tion of the continuum to 25 fim is in reasonable agreement 
with the observed flux density; however the extrapolation 
to 5 fim fails to predict the 4-5 fim photometry of IRAS 
05189—2524 (see double arrow), suggesting that an unre- 
lated power source is responsible for this emission. Given 
the apparent similarity in the shape of 1-5 fim emission in 
both IRAS 05189-2524 and IRAS 13349+2438, it is possible 
that the spectral energy distribution of IRAS 05189—2524 is 
the result of a starburst producing the majority of the 8-100 
fim emission and an AGN with a spectral energy distribu- 
tion similar to that of IRAS 13349+2438 responsible for the 
bulk of the shorter wavelength emission but contributing lit- 
tle to the 8-100 fim emission. Such a decomposition could 
reconcile the apparently contradictory observations of both 
an AGN-like optical spectrum and the present PAH type 
8-13 fim spectrum. A decomposition of this kind could also 
be consistent with the suggestion of Goldader et al. (1995) 
that the 2 fim continuum is dominated by AGN [heated 
dust] emission based on the small equivalent width of the ob- 
served Br7 emission. That hot dust typically plays a lesser 
role in starburst galaxies at 2 fim may also be deduced from 
the frequent detection of stellar continuum with probably 
little veiling of the CO absorption features, while in IRAS 
05189—2524 CO absorption is not detected (Goldader et al. 
1995). Another piece of evidence that supports the view that 
the AGN in IRAS 05189-2524 may not be energetically 
dominant is the failure to detect a radio core of extreme 
brightness temperature in VLBI observations (Lonsdale, et 
al. 1993). 

The VLA radio map of CHYT shows the radio con- 
tinuum in this source to be extended on a scale of 160 pc, 
which is comparable to the minimum blackbody diameter of 
280 pc based on the observed 60 and 100 fim flux densities 
and assuming optically thick dust emission. It is therefore 
plausible that the far-infrared emission arises from optically 
thick warm dense clouds embedded in the starburst with a 
low volume filling factor but a column covering factor of or- 
der unity. The mid-infrared continuum and PAH emission 
present in this spectrum might arise from the surfaces of 
these clouds and may be optically thin in the sense that the 
emission arises from cloud surfaces facing our line of sight, 
and from a layer that is only 1-2 Ay thick on any given line 
of sight into the starburst. 

Given the possibility of a decomposition of the 1- 
100 fim spectral energy distribution into a dominant star- 
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burst component and an AGN component important only at 
shorter wavelengths, the plausibility of the infrared emission 
being spatially associated with the radio continuum, and 
most importantly, the presence of 11.3 /im PAH emission in 
the 8-13 /im spectrum taken together with the adequate fit 
of the model given in Fig. 1, the 8-1000 pm emission from 
IRAS 05189-2524 is attributed to a starburst for the pur- 
poses of this work. 

4-2.2 Digression on newly available ISO data 

The model applied to the data presented in Fig. 1 for the 
purpose of spectral classification is not required in cases 
where the strength of the 7.7 fim PAH feature has been 
observed as is now the case for four of the ultraluminous 
infrared galaxies in this combined sample as a result of In- 
frared Space Observatory (ISO) observations (Genzel et al. 
1998). With the additional information it is possible to com- 
pare observations with a variety of scenarios to better un- 
derstand the physical situation with greater confidence than 
possible when only ground-based data are available. In Fig. 
3, three possibilities are considered: (1) For the top spectrum 
in each panel, a PAH plus power-law continuum model such 
as that employed in Fig. 1 is fit not just in the spectral 
region surrounding the 11.3 fim feature, but to the entire 
6-13 fim range. (2) For the middle spectrum in each panel, 
a PAH plus power-law continuum model is subjected to ex- 
tinction by cold dust that includes the silicate absorption 
feature, as adopted by Dudley & Wynn- Williams (1997), 
and again fit to the entire 6-13 /jm range. (3) For the bottom 
spectrum in each panel, a PAH plus power-law continuum 
model is additively combined with a silicate-absorbed con- 
tinuum source, and again fit over the entire 6-13 fim range. 
The model fits presented in Fig. 3 have been carried out 
using CURVEFIT with equal weighting. In the first model, 
the index of the power-law was specified either from Table 
3 (IRAS 17208-0014 and Mkn 273) or as 12.8 and 4.9 for 
Arp 220 and Mkn 231, respectively. It should be noted that 
three of the sources exhibit a large dynamic range, and in 
the following discussion some attention will be drawn to the 
poor fit at low flux density levels, which is also an expected 
systematic effect of the fitting procedure, since differences 
between small numbers are small. It may also be noted that 
the mismatches are apparent for a large number of neigh- 
boring data points. 

4-2.3 IRAS 17208-0014 

In IRAS 17208-0014, all models give fairly adequate fits 
to the data, and deciding between them depends mainly on 
other observational evidence. The strongest evidence that 
supports a preference for the fit shown for the top spectrum 
is that this source is extended on a scale ~1.7 kpc in the ra- 
dio (Condon et al. 1996), although higher spatial resolution 
observations reveal the presence of a sub-kpc radio contin- 
uum source associated with the OH maser emission (Martin 
et al. 1989). The failure of Martin et al. (1989) to detect 
more extended emission may be due to the restricted set 
of base lines they employed. When this evidence is taken 
along with the optical spectral classification of Veilleux et 



al. (1995) in which the line ratios in this source are similar 
to those in Hll regions, there is no strong reason to suspect 
that anything other than star formation is responsible for 
the observed mid-infrared spectrum. Nor is there any rea- 
son to suspect a high degree of extinction; IRAS 17208—0014 
is among the sources presented by Genzel et al. (1998) for 
which the continuum underlying the 12 A- fim [Ne II] feature 
in their SWS data can be comfotably plotted on the scale 
of the ISOPHOT-S data in their fig. 2, suggesting a lack of 
strong silicate absorption (see also Fig. 1). The fit to the 
upper spectrum might be improved somewhat by increasing 
the strength of the 7.7 and 6.2 /jm features by a factor of 2, a 
procedure which is probably allowed for the 7.7 fim feature 
at least since the Orion Bar spectrum used as a template 
here is known to have a somewhat weaker than usual 7.7 
(im feature relative to the 11.3 fim feature (Zavagno et al. 
1992). 

4.2.4 Mkn 273 

The fit to the top spectrum of Fig. 3 for Mkn 273 gives, as 
might have been anticipated from Fig. 1, a poor representa- 
tion of the data. The fitted strength of the 11.3 (im feature 
is clearly too large, which is roughly the converse of the re- 
sult in Fig. 1. Thus, it may be concluded that a simple PAH 
plus power-law continuum model does not describe these 
data. The fit shown for the middle spectrum is better, in 
that it requires less PAH emission over all and the strength 
of the 11.3 fim feature is not over-predicted; however the 
continuum around the 9.7 fim silicate minimum is not well 
reproduced. This suggests that an unaccounted-for emission 
component is required to produce this emission. The fit to 
the lower spectrum overcomes this difficulty by superimpos- 
ing an unobscured PAH plus power-law continuum such as 
that which fit IRAS 17208-0014 (top spectrum) upon a con- 
tinuum source that suffers strong silicate absorption. Since 
Mkn 273 shows clear signs of nuclear activity (Armus, Heck- 
man & Miley 1989; Veilleux et al. 1995), a silicate-absorbed 
continuum source might be produced by an edge-on view of 
a disk model such as those proposed by Pier & Krolik (1992) 
or Efstathiou & Rowan- Robinson (1995), which are natural 
elaborations of the AGN unification model proposed by An- 
tonucci & Miller (1985). The continuum level evident in the 
SWS observations of 12.8-fim [Nell] in Genzel et al. (1998) 
is consistent with the new CGS3 data presented here; and 
therefore consistent with silicate absorption playing a role 
in this source. If the fit to the bottom spectrum is accepted, 
then the AGN component could well be responsible for the 
greater fraction of the 60 and 100 fim emission. 

4.2.5 Arp 220 

Prior to the new ISO observations the analysis of the Arp 
220 spectrum would have proceeded very much along the 
lines that have just been given for Mkn 273 (Smith et al. 
1989), with the proviso that Arp 220 has less obvious signs 
of activity in its nucleus so that disk-like models for the 
silicate-absorbed continuum emission may be less helpful in 
advancing our understanding, and deeply embedded models 
perhaps more descriptive (Dudley & Wynn- Williams 1997). 
Applying the model used in Fig. 1 to the data of Smith et 
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al. (1989) severely under-predicts the observed 8.2 /im flux 
density. However, as can be seen in Fig. 3 the new ISO data 
can apparently be fit by a simple PAH plus power-law con- 
tinuum model due to the apparent discrepancy between the 
observations of the 11.3 /im feature reported by Smith et al. 
(1989) and the ISO data of Genzel et al. (1998). It should 
be noted that there are no continuum points longward of 
the 11.3 fim feature in the new ISO data, so that determin- 
ing its strength from those data is difficult. It seems safer 
at this point to accept the ground-based data as a guide in 
this spectral region and thus hold this model in abeyance, 
pending further observational evidence. The fit to the mid- 
dle spectrum is essentially the same as that given in fig. 1 
of Smith et al. (1989) and suffers from the same difficulty 
noted by Smith et al. (1989), namely that the minimum 
near 9.7 fim is poorly fit. The fit to the lower plot is similar 
to that shown in Smith et al.'s (1989) fig. 2 or as adopted 
by Dudley & Wynn- Williams (1997). The continuum in the 
12.8- fim [Nell] SWS data of Genzel et al. (1998) is consis- 
tent with the Smith et al. (1989) data but is not apparently 
consistent with the ISOCAM-CVF spectrum presented by 
Elbaz et al. (1998). 

For Arp 220, therefore, the new ISO data presented by 
Genzel et al. (1998) seems to allow the possibility that unex- 
tinguished PAH emission plus a simple power-law continuum 
describes this spectrum as a result of the apparent strength 
of the 11.3 fim feature in the new data. It has been suggested 
that differing spectral resolution (roughly a factor of 2) be- 
tween the Smith et al. (1989) and Genzel et al. (1998) data 
may account for this difference, however this could only work 
and if the 11.3 (im feature were unresolved in both spectra, 
and that does not appear to be the case. Gauging the over 
all PAH emission from the 6.2 fim feature which displays 
continuum on either side, the intensity of which is found to 
be well correlated with the that of the 7.7 /jm feature in 
Galactic sources (Cohen et al. 1989), it is clear the PAH 
emission relative to the over-all 8-1000 fim emission is ex- 
ceedingly weak, really no stronger than what might be al- 
lowed in Mkn 231 (see next), if the apparent 6.2 /jm feature 
in the ISO data for Mkn 231 is real. In contrast, again using 
the 6.2 /jm PAH feature as a guide, the strength of the PAH 
emission relative to the 8-1000 /im emission in Mkn 273 
is twice as strong, and in IRAS 17208—0014, five times as 
strong. The Arp 220 spectrum might possibly be produced 
by a starburst that reprocesses more that 90 per cent of its 
luminosity not once but many times so that the observed 
PAH emission is produced only in a thin outer shell or in 
a small number of narrow chinks. Such a model would be 
even more extreme than what is argued above as plausible 
for IRAS 05189—2524 in a starburst context. However, given 
the disagreement between the two spectra of this source near 
11.3 /im and in consideration of the arguments of Dudley & 
Wynn- Williams (1997), a model that has both unobscured 
PAH emission and a silicate-absorbed continuum due to a 
deeply embedded AGN continues to seem more natural. This 
interpretation also seems to account simultaneously for both 
the 158-fim [CII] (Luhman et al. 1998) and 2-10 keV X-ray 
(Iwasawa 1999) deficits in this source if the "AGNs" (Soifer 
et al. 1999) are sufficiently obscured. It should be noted that 
in this preferred model, only a fraction of the flux density 
at 7.7 fim can be attributed to the 7.7 fim PAH feature, 
in accordance with either the 12.8 of 14.3 /mi continuum 



seen in the SWS data presented by Genzel et al. (1998). 
As suggested by Dudley & Wynn- Williams (1997), spatially 
resolved spectroscopic observations may help to clarify this 
issue. 



4.2.6 Mkn 231 

For Mkn 231 any PAH contribution to the spectrum must 
be minimal compared to a continuum contribution (Roche 
et al. 1983), as can also be seen from the fit to the upper 
spectrum in Fig. 3. This source had been fit adequately with 
a model consisting of a continuum source suffering silicate 
absorption by Roche et al. (1983), and the present models 
are not an improvement. It would be of interest to know 
if disk radiative transfer models such as those proposed by 
Pier and Krolik (1992) would fair as well as the cold ab- 
sorber model of Roche et al. (1983) in the detailed fit to the 
silicate absorption feature given their likely success in fitting 
the overall infrared spectral energy distribution of Mkn 231. 
The 11.3 fim emission feature is not detected in this source 
and it is not clear that the single point that is spectrally co- 
incident with the 6.2 fim feature in the ISO data should be 
considered as a real indication of the presence of PAH emis- 
sion. Some PAH emission is not unexpected in this source 
given the presence of star-forming knots near the nucleus of 
this source that might account for ~10 per cent of its in- 
frared luminosity (Surace et al. 1998). The estimate of the 
7.7 fim PAH strength given by Genzel et al. (1998) relies 
on a continuum point at 10.9 fim which lies in the silicate 
absorption feature. Of the sources measured in this way by 
Genzel et al. (1998) Mkn 231 perhaps provides the clear- 
est example of how this method is prone to overestimate 
the PAH strength when silicate absorption is present. The 
12.8 or 14.3 fim continuum in their fig. 1 provide a rough 
guide as to when there may be a problem in this respect. 
For Mkn 231, the 12.8 fim continuum of Genzel et al. (1998) 
is high in comparison with the data of Roche et al. (1983) 
but it is possible that it would agree with the extrapolation 
of the long-wavelength ISOPHOT-S data which also appar- 
ently fall significantly above the Roche et al. (1983) data as 
seen in Fig. 3. It is not clear if these differences should be at- 
tributed to aperture effects, or other systematics, or if they 
reflect intrinsic source variability over 1.4 decade; Roche et 
al. (1983) found good agreement with the narrowband pho- 
tometric observations of Reike (1976) made over 0.4 decade 
earlier. 



4.3 A 60 fim flux-limited subsample 

Among galaxies with Lir> 1.6 x 10 11 Lq, a nearly complete 
(all but three sources have 8-13 fim spectra) flux-limited 
(at 60 fim) sample of 25 sources may now be defined on the 
combined Bright Galaxy Survey (Soifer et al. 1989) and the 
Extended Bright Galaxy Survey (Sanders et al. 1995) (com- 
bined BGS) for F„(60) > 11.25 Jy (roughly twice the flux 
limit of combined BGS). The only other restriction on this 
subsample is that the sources be observable using UKIRT, 
i.e. —25° < S < 60°. Table 5 gives the ultraluminous and 
luminous samples with the source name in column 1, the 
8-13 spectral classification in column 2, the reference for 
the 8-13 fim spectrum in column 3, the log of the 8-1000 
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/mi luminosity in column 4, and the log of the ratio of the 12 
fim flux density to the 60 /jm flux density in column 5. The 
ultraluminous sample has 5 members and the luminous com- 
parison sample has 20 members. The ultraluminous sample 
is unlikely to span all the possible characteristics of ultralu- 
minous infrared galaxies and therefore is not fully represen- 
tative. There are two observations that confirm this: (1) if 
Arp 220 were scaled to the 60 /jm 11.25 Jy flux limit used to 
define this sample, it would not have been observable using 
CGS3 or other ground-based spectrographs available when 
these data were collected, and (2) the ratio given in column 5 
of Table 5 for the two brightest ultraluminous galaxies (Arp 
220 and Mkn 231) bracket the same ratio for all the galax- 
ies in the sample with 1.6 x 10 11 < Lir< 10 12 Lq except 
NGC 1068 and NGC 7469. An analysis over the entire Bright 
Galaxy Survey shows that the dispersion of this ratio is in- 
deed larger for ultraluminous infrared galaxies than for lumi- 
nous infrared galaxies by a factor of 2 (Soifer & Neugebauer 
1991). Thus, conclusions based on the present flux- limited 
sample of ultraluminous galaxies should be viewed with cau- 
tion not simply because of the small number of galaxies in 
this sample (counting statistics on a bimodal parent popu- 
lation), but because the characteristics of the parent popu- 
lation are not well represented (under-sampling). 

It is possible that the last of these two difficulties could 
be somewhat ameliorated by the following considerations. 
Ultraluminous galaxies that fall outside the range of 12 on 
60 /im flux density ratios defined by Arp 220 and Mkn 231 
may be assigned to the AGN-dominated bin in the context 
of exploring the AGN vs. starburst hypotheses. For ultralu- 
minous galaxies with 12/60 /im colors redder than Arp 220, 
PAH emission cannot make a contribution to the 12 /im flux 
density that is sufficiently large to explain the 60 /im flux 
density. It is therefore likely that such sources will have 
deep silicate absorption features, and by the arguments pre- 
sented by Dudley & Wynn- Williams (1997), the size of their 
power source should be too small to be due to a starburst. 
For galaxies with 12/60 /im colors bluer than Mkn 231, 60 
/im selected samples will overlap 25 or 12 /im selected sam- 
ples that are known to select for AGNs (de Grijp et al. 1985; 
Spinoglio & Malkan 1989). However, reliance on these no- 
tions would be ill advised given the internal evidence in the 
present data showing that PAH emission can be important 
in ultraluminous galaxies with 12 to 60 /im flux density ra- 
tios similar to the two extremes, namely, IRAS 17208—0014 
and IRAS 05189-2524. 

Inspection of Table 5 shows that 3 out of 5 of the ul- 
traluminous galaxies are powered predominantly by AGNs 
although if judgment is reserved for Arp 220, this fraction 
would be 2 out of 5. The trend is thus consistent with about 
half of ultraluminous galaxies being powered by AGNs. On 
the other hand, the fraction of the more numerous luminous 
galaxies in this sample that are powered predominantly by 
AGNs is at most 0.25 and could be as small as 0.05 or smaller 
if the relative contributions of the AGN and starburst in 
NGC 1068 are approximately equal (Telesco et al. 1984). 
This trend supports the evidence from optical spectroscopy 
that AGNs are present in greater numbers in ultraluminous 
infrared galaxies (Veilleux et al. 1995). It goes further than 
optical or radio studies because it is possible to demonstrate 
with greater certainty which type of power source (starburst 
or AGN) dominates the infrared luminosity. Further, it is 



apparent that for luminosities less than 10 Lq, star for- 
mation is far and away the dominant power source, a re- 
sult that confirms work at many wavelengths. This analysis 
lends some support to the proposal of Sanders et al. (1988) 
that there is an evolutionary connection between ultralumi- 
nous infrared galaxies and quasars, while at the same time 
demonstrating that massive star formation can indeed be 
the dominant power source of some ultraluminous infrared 
galaxies. However, it must be conceded that the evidence 
presented constitutes only a trend. The precise proportion 
of AGN-dominated sources among 60 /jm selected ultralu- 
minous infrared galaxies cannot be said to have been estab- 
lished here, and thus no statistically meaningful difference 
between luminous and ultraluminous galaxies has truly been 
shown. It is of interest that the trend does continue with 
the next step down in 60 /im flux density. Adding IRAS 
12112+0305 and IRAS 08572+3915 to the present sample of 
five ultraluminous galaxies adds one apparently starburst- 
dominated (based on the apparently normal ratio of the 6.2 
/xm PAH feature emission to overall infrared emission in the 
IRAS 12112+0305 spectrum of Genzel et al. 19980) and one 
AGN-dominated object (Dudley & Wynn- Williams 1997). 



4.4 Anti-correlation between 11.3 /im eqivalent 
width and L m /M(R 2 ) 

In Fig. 4 a possible anti-correlation between the equivalent 
width (EW) of the 11.3 /im feature and the global ratio of 
Lm/ M{R2) is presented for a number of galaxies. The EW 
data are based either on the model fits given in Table 3 
or the previously published values of RASW reproduced in 
Table 6. The ratios of infrared luminosity to H2 mass are 
taken from the literature (Tables 4 and 6). The solid line in 
Fig. 4 indicates inverse proportionality. 

The ratio , is often taken to be a measure of star 
formation efficiency. When it is large, then a large fraction 
of the available gas is involved in star formation; when it is 
small, conditions may be such that the gas exists in a more 
quiescent state. However, for a very large ratio of ; , 
Sanders et al. (1991) argued that even a very high star for- 
mation efficiency cannot account for the infrared luminosity, 
on the grounds that star formation would too quickly con- 
sume the available gas mass, so that AGN-like activity might 
be the more likely primary power source of the infrared lu- 
minosity. Both ; and EW(11.3) are independent of dis- 
tance in so far as EW(11.3) is independent of aperture size. 

A tendency for EW(11.3) to be anti-correlated with 
the star formation efficiency might arise if two conditions 
were met. First, if the intensity of the 11.3 /im feature is 
a roughly constant fraction of the overall infrared emission 
(proportional to Lm) then this will not contribute to the 
anti-correlation. Second, if the strength of the continuum at 
11.3 /im increases with decreasing availability of dust (in- 
versely proportional to M(H2)) to reprocess the light emit- 
ted by massive stars so that a larger fraction of the dust 

t A crude attempt to correct for the contribution of the silicate 
absorption edge to the 7.7 fim PAH feature by fitting a power law 
to continuum points at 5.9 and 14.3 /im suggests a deficit in the 
strength of the 7.7 /im feature for IRAS 12112+0305. 



© 1999 RAS, MNRAS 000, 000-000 



Infrared galaxies 9 



grain population is heated to the warm temperatures re- 
quired for the dust grains to emit at 11.3 /im then this could 
be the source of anti-correlation. Such an explanation re- 
quires further elaboration since the process by which the 
11.3 /im continuum is produced is likely, in general, to be a 
local phenomenon rather than a global one. 

An elaboration that may be in keeping with recent HST 
observations (Watson et al. 1996; Surace et al. 1998; Satya- 
pal et al. 1997) is that as the star formation efficiency in- 
creases, the proportion of massive stars formed in dense clus- 
ters may increase. If this is the case, then the time-scale 
for clearing the interstellar medium from around these clus- 
ters to a distance where dust grains would no longer be 
heated sufficiently to contribute to the 11.3 /jm continuum 
may be longer that the time-scale for clearing the interstellar 
medium around the groups of massive stars that might form 
if the star formation efficiency is lower. Two effects may con- 
tribute to this. First, the distance at which a grain could be 
heated to a temperature high enough to contribute to the 
11.3 /im continuum would increase as the square root of the 
number of massive stars in a cluster so that if the velocity 
of gas and dust being cleared from clusters of massive stars 
is not too different from the velocity of gas and dust cleared 
from groups of massive stars, the dust remains heated for a 
longer period of time. Second, dense clusters might be tem- 
porarily gravitationally bound and behave ballistically with 
respect to the dissipative gas. They may sometimes move 
into a fresh reservoir of gas, heating the associated dust to 
higher temperatures. 

In this view, the additional contribution to the 11.3 
/tm continuum emission at higher star formation efficiency 
would be provided by possibly equilibrium temperature 
grains heated by a combination of direct and Lya heat- 
ing rather than small grains heated by a single photon. To 
zeroth order, this is a sensible proposition, since one ex- 
pects (at least within photo-dissociation regions) that the 
PAH emission and the spike-heated grain emission to have 
roughly fixed relative contributions at 11.3 /im due to the 
rough similarity of the emission time scales for the two emis- 
sion processes, so that the apparent anti-correlation would 
not necessarily arise. 

A number of observations will be required to examine 
this situation. Most important are higher signal-to-noise ob- 
servations that will allow real measurements of EW(11.3) 
for a number of the sources included in Fig. 4. These are 
required to decide if the anti-correlation is has any merit. 
Another avenue of investigation would be to make higher 
spatial resolution observations of nearby starbursts where 
clusters of massive stars are observed to exist to estab- 
lish whether or not EW(11.3) does decrease in the vicin- 
ity of some clusters (see Dudley & Wynn- Williams 1999 for 
a preliminary attempt). Finally, high spatial resolution 2 
jjm adaptive optics imaging of a number of starbursts with 
a range of star formation efficiency can help to determine 
if the fraction of massive stars that exist in dense clusters 
varies with the star formation efficency parameter as con- 
jectured here. 



5 CONCLUSIONS 

New 8-13 fim spectroscopy of 27 infrared galaxies has been 
presented. The spectra have been fit with a two-component 
model consisting of a power-law continuum underlying a 
PAH emission spectrum. These fits are based on the spec- 
tral data between 10 and 13 /jm and are used to predict the 
flux emitted between 8.0 and 8.4 /im. It is found that for 
the majority of the observed galaxies, the model predicts 
the observed 8.0-8.4 /xm emission to within an acceptable 
factor of 2, suggesting that the majority of the sources are 
powered by starbursts. 

One source, Mkn 273, has been found to show very clear 
evidence for silicate absorption. 

A flux-limited 60 /mi selected subsample that combines 
the present spectroscopic results and those of RASW had 
been defined, consisting of 25 galaxies. In this subsample 
the fraction of ultraluminous infrared galaxies powered by 
AGN is ~50 per cent in contrast to 5-25% among the lumi- 
nous infrared galaxies. Owing to the small sample size for 
the ultraluminous infrared galaxies, this difference in AGN 
fraction can be considered only a trend. 

A possible anti-correlation is identified between the 
equivalent width of the 11.3 fim PAH feature and the star 
formation efficiency in galaxies where both have been mea- 
sured. It is suggested that variation in the continuum dust 
temperature as a function of star formation efficiency could 
produce this apparent effect. 
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Figure 1. 8—13 /im spectra of 27 galaxies are presented. The small filled circles are the 
(nearly) independent data points discussed in Section 2.1 with 1 a error bars. The filled 
circles, which are either twice as large or four times as large as the spectral data, are the 
same data binned over the ranges indicated by the horizontal bars. The largest of these 
filled circles are shown only when one or both of the rest wavelength bins, 10—10.5 or 
10.5—10.9 /im, represent detections of <1 a significance. Apparent differences between the 
spectral and binned data are due to missing spectral data points that, by themselves, are 
upper limits only. The dashed lines represent the best-fitting power-law long- wavelength 
continuum underling a PAH spectrum of the Orion Bar represented by solid or dot- 
dashed lines. Dot dashed lines are used when the formal uncertainty in the strength of 
the fitted PAH component is more than three times greater than the fitted strength of the 
component and are shown at the formal 3 a upper limit to the strength PAH component. 
Diamonds represent IRAS 12 /im fluxes, the sources of which are given in Table 4. For 
Arp 299B1, the high-resolution spectrum is also presented in the figure divided by a factor 
of 10 for clarity. 

Figure 2. The spectral energy distribution of IRAS 05189 — 2524 is compared with those 
of another ultraluminous infrared galaxy: Mkn 231 and the infrared loud quasar IRAS 
13349+2438. In the figure, 0.9—5 /im photometry (open diamonds) are from McAlary, 
McLaren & Crabtree (1979) for Mkn 321; from Dudley (1998) for IRAS 05189-2524 with 
the exception of the 0.9 /im point which, is from Surace et al. (1998); and from Beichman 
et al. (1986) for IRAS 13349+2438 (including the IRAS data). 8-13 spectroscopy 
(filled circles) is from Roche, Aitken &; Whitmore (1983) for Mkn 231 and from this 
work for the remaining two sources. IRAS data are given by open squares. The data for 
IRAS 13349+2438 have been presented in the quasar's rest frame, and the spectral energy 
distributions of IRAS 05189-2524 and IRAS 13349+2438 have been divided by factors 
of 2 and 20, respectively. 



Figure 3. The combined ground-based and ISO spectra of four ultraluminous galaxies 
are presented in the observed frame. For each galaxy, the top spectrum gives the observed 
spectrum, while the lower two are shifted by multiplicative factors. The ground based 
data (filled circles) are from this work for IRAS 17208-0014 and Mkn 273, from Smith 
et al. (1989) for Arp 220, and from Roche et al. (1983) for Mkn 231. The ISO data (open 
circles) are from Genzel et al. (1998) and have been corrected by a factor of 1.4. Three 
dust emission models are shown for each source with the total emission given by the solid 
lines and continuum emission given by the broken lines; for the fits to the bottom spectra, 
an additional unobscured power-law continuum component is not shown. 

Figure 4. A possible anti-correlation between star formation efficiency and the equivalent 
width of the 11.3 /im PAH feature is shown. 
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Table 1. Observing log and adopted redshifts for spectroscopic observations. 
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Table 1 — Continued 



Object 


a(1950) 


6(1950) 








Date 




App 


Chop 


PA 


Name 


(h m s) 


r ' ") 


Ref 
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Ref 


(UT) 
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1995 Aug 
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IRAS 20210+1121 


20 21 02.58 


+ 11 21 52.5 
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0.0564 
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1994 Aug 
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15 


135 














1995 Aug 
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II Zw 96 


20 55 05.3 


+ 16 56 03 
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0.0363 
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1994 Aug 


29 


3.88 
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45 














1995 Aug 


14 


3.26 


20 


90 


Zw 475.056 


23 13 33.13 


+25 17 01.6 
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0.0274 


1 


1993 Aug 


31 


2.40 


15 


90 














1993 Sep 01 


3.26 







References. — Positions are from (1) this work, (2) CHYT, (3) Hewitt & Burbidge (1989), (4) Eales et al. 
(1990), (5) Condon et al. (1996), (6) Wynn- Williams et al. (1991), (7) Condon et al. (1990), (8) Condon, Frayer 
& Broderick (1991), and (9) Clements (1981). Redshifts are from (1) Soifer et al. (1989), (2) Strauss et al. (1992), 
and (3) Sanders et al. (1995) 

a Brightest optical component 

b Brightest K-band component 

c From figure rather than table 

d 8-13 /tin peak up position 

e The redshift for this source has been dropped from later versions of this catalog. Crawford et al. (1996) failed 
to detect this source at radio wavelengths, and propose that it may be attributable to Galactic infrared cirrus. The 
present spectrum (see Fig. 1) is not typical of cirrus, so the source has been retained here. The reported sensitivites 
of Crawford et al. (1996) suggest that this source is under-luminous at 20 cm by a factor of at least 8 compared 
to the usual far infrared-radio correlation (Helou, Soifer, & Rowan- Robinson 1985). 



Table 2. Calibration stars 



Object 


Alternate 


Date 


HR 


F„(10.1) 




T 


Name 


Name(s) 


(UT) 


(#) 


(Jy) 


Ref 


(K) 


Mkn 551 


IV Zw 023 


1995 Aug 14 


7001 


39.8 


1 


11000 


IC 1632 


Arp 236, ESO 541-IG 023 


1992 Oct 16 


8775 


413 


1 


3000 


NGC 695 


UGC 01315, V Zw 123 


1993 Aug 30 


7001 


39.8 


1 


11000 


NGC 828 


UGC 01655, VI Zw 177 


1992 Oct 14 


2491 


147 


1 


10000 






1992 Oct 15 


8775 


413 


1 


3000 


MCG +05-06-036 


V Zw 233, Mkn 1034 


1994 Aug 29 


0188 


57.5 


3 


5800 


IRAS 02528 + 4350 




1993 Aug 31 
1993 Sep 01 


7001 


39.8 


1 


11000 


IRAS 04154 + 1755 




1995 Mar 19 
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238 


1 


5500 






1995 Mar 20 


1457 


648 


1 


5000 






1995 Mar 21 


1457 


648 


1 


5000 


IRAS 05189-2524 




1993 Feb 11 


2943 


80.1 


1 


8000 


NGC 2388 


UGC 03870 


1992 Feb 16 


3748 


150 


1 


4000 


NGC 2623 
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1992 Feb 17 


4518 


18.1 


2 


6000 
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1995 Mar 19 
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4000 
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15.4 


2 


5000 


Arp 299C a 
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4518 


18.1 


2 


6000 


Arp 299B1 


NGC 3690 


1991 May 17 


5340 


738 


1 


4400 


Arp 299C' 




1995 Mar 20 


4518 


18.1 
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6000 


Arp 299A 


IC 0694 


1991 May 18 


4534 


7.31 


1 


8700 
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4518 


18.1 
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6000 
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5235 


11.5 
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6000 
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1 


4400 


IRAS 17208—0014 




1993 Aug 30 
1993 Aug 31 

1993 Sep 01 

1994 Aug 29 


7001 


39.8 


1 


11000 


IRAS 18293-3413 




1993 Sep 01 


7001 


39.8 


1 


11000 


ESO 593-IG008 




1994 Aug 27 


7264 


8.32 


3 


9000 






1995 Mar 20 


1457 


648 


1 


5000 






1995 Aug 13 


7557 


33.3 


2 


9000 


IRAS 20210 + 1121 




1994 Aug 28 


7949 


39.8 


3 


6000 






1995 Aug 13 


7557 


33.3 


2 


9000 


II Zw 96 




1994 Aug 29 


7949 


39.8 


3 


6000 






1995 Aug 14 


7001 


39.8 


1 


11000 


Zw 475.056 


IC 5298, MCG +04-54-038 


1993 Aug 31 
1993 Sep 01 


7001 


39.8 


1 


11000 



a This spectrum was previously reported in Dudley Sz Wynn- Williams (1993); it has now been recalibrated to be ~10% fainter. 



References. — (l) IRTF Photometry Manual, (2) J. K. Davies (1996; private communication), (3) assumed 



Table 3. Spectral measurements of galaxies. 



Object 


a 


1(11.3) 


EW(11.3) 


/„ S ource(8.2) 
Sv model(8.2) 


I([SIV]) 


I([Nell]) 


Name 




(1CT 15 W m" 2 ) 


(/tin) 


(1CT 15 W m" 2 ) 


(1CT 15 W m" 2 ) 


TV/Tim r ci 

IVLkn 551 


ii a _i_n a 
ll.bzbU.b 


<2 
0.85+0.06 

1.610.3 
1.56+0.04 

0.810.2 
0.54+0.03 


0.8 + 1.9 


i a _l_n o 


<u.y 




1C Iboz 


o o in o 
o.ozbU.Z 


0.43+0.25 


2.2±U.2 


A A 

<u.y 


<1 


at a o k 
IN by 5 


y .ozbu.y 


1.4+4.5 


n o _l_n o 
U.o±U.2 


<U.O 




AT f" 1 f" 1 Q O Q 
IN KjLj ozo 


o a I n o 
Z.4±U.O 


<2 
0.84+0.05 


1+2 


1.2±U.3 


<1 


<z 


MCG +05-06-036 


1.4+0.6 


<1 
0.77^0.05 

<0.4 


0.9+2.0 


0.8±0.2 


<0.9 




ll\Ao Uz5zo+4o5U 


n £ _i_n o 
U.bzbU.o 


-0.012+0.023 


<u.y 


<5 




ll\Ao 1)4154+1(55 


5.o±U.o 


<0.9 
0.55+0.04 


0.34+0.33 


i i _i_n o 
l.lzbU.Z 


A r 
<U.5 




IRAS 05189-2524 


2.47+0.05 


<2 
2.05+0.05 
2.5+0.3 

3.11+0.04 


0.20+0.04 
0.6 + 0.1 

0.94+0.49 


1.02 + 0.09 


<1 




IN IjC zooo 


£ n_i_n i 
b.UzbU.l 


U.o5±U.U ( 


<U. ( 


i o_i_n o 
l.ZzbU.o 


IN IjC ZD ZO 


in c±n a 
lU.b±U.4 


<0.9 
0.76+0.04 


0.52+0.70 


o i _Ln o 
Z.lzbU.o 


<U.b 


n o i n o 
U.ozbU.o 


ttj a o nm 1 1 i nn7 
ll\Ao Uylll-lUU( 


o n_Ln a 
o.U±U.4 


<0.6 
0.16+0.03 


0.17 + 0.27 


l.b±U.5 


<U.4 




TTJ A O nO/l 1 C01 

ll\Ao Uy4oo-15ol 


( .o±U.5 


<0.7 

0.54+0.03 
1.5+0.3 

1.46+0.04 


0.59+0.98 


i o_i_n o 
l.ZzbU.Z 


A A 
<U.4 




at 01 1 n 
IN ollU 


K i _Ln c 
5.1±U.5 


1.5+2.5 


n o i n i 
U.ozbU.l 


<U. 1 


A Q 
<U.O 


TTJ A O 1 nC^C 1 O A A Q 

ll\Ao lU5b5 + z44o 


ii i _Ln o 
ll.lzbU.o 


<0.7 
0.73+0.04 


0.40+0.41 


o a _Ln o 
Z.4±U.Z 


<U. ( 




Arp zyyc 


4.zU±U.Ub 


2.0+0.2 
1.99+0.04 


0.27 + 0.04 
0.28+0.06 


l.U ( zbU.Ub 


<U.b 


o n_Ln o 
Z.UzbU.Z 


Arp zyyjDi 


b.b4±U.U5 


<2 
1.21+0.06 


0.079+0.014 


o o _i_n o 
o.ozbU.Z 


A r a 

<U.5 


n 7±n o a 
U. ( ±U.z 


Arp 299C 


9.67+0.01 


<0.9 
0.54+0.04 

4 + 1 
5.10+0.08 
1.4+0.5 

2.29+0.05 


0.27 + 0.02 


2.7+0.3 


<0.7 


<0.9 


Arp zyy 


o o i _l_n no 
o.21±U.U3 


0.13+0.04 
0.16+0.02 


o o _l_n i 


<2 


^- A 

<4 


Arp zyyA 


ik o _l_n i 
15.o±U.l 


0.41+0.19 


2.D±U.1 


^ n o 
<U.o 


<z 


T T / 1 / 1 AOOQ7 


ill _l_n /i 
11.1±U.4 


0.9+0.2 
0.65+0.03 


0.62+0.93 


o a _l_n o 
2.4±U.2 


^ d A 

<U.4 


n n_Ln o 

u.y±u.z 


TTJ A C 100/1O I O/IOO 

ll\Ao loo4y + z4oo 


i o o _l_n no 


<0.8 
<0.1 


-0.020+0.008 


i nc i n n a 
l.U5±U.U4 


^ n k 
<U.5 




T T /" 1 /" 1 (\oac\G. 

U IjC Uobyb 


in /i _l_n o 
1U.4±U.2 


<0.8 
0.35+0.04 


0.16+0.14 


a n_Ln a 
D.U±U.D 


^ n v 
<U. 1 




I Zw 107 


6.7±0.3 


<0.7 
0.60+0.04 


0.38+0.36 


1.6±0.2 


<0.5 




NGC 6240 


9.8±0.2 


2.9+0.6 

2.20+0.05 
1.2+0.3 

1.54+0.04 


0.47 + 0.26 


1.4±0.2 


<1 


<2 


IRAS 17138-1017 


8.7±0.2 


0.56+0.40 


1.4±0.1 


<0.8 


<0.8 


IRAS 17208-0014 


7.8±0.2 


<1 
0.77 + 0.04 


0.40+0.35 
0.8 + 0.3 

0.83+0.43 


1.7±0.1 


<0.6 




IRAS 18293-3413 


4.4±0.1 


3.7 + 0.5 
3.35+0.05 


0.7±0.2 


<1 


<1 


ESO 593-IG008 


6.7±0.6 


<1 
0.47 + 0.05 


0.7 + 1.4 


1.4±0.3 


<0.6 




IRAS 20210+1121 


3.2±0.1 


<1 
0.58+0.05 


0.11+0.04 


0.8±0.1 


<0.7 




II Zw 96 


6.9±0.2 


<1 
<0.1 


-0.015+0.022 


2.4±0.5 


<0.7 




Zw 475.056 


2.0±0.1 


<0.9 
0.82+0.04 


0.25+0.13 


1.0±0.1 


<0.7 


<1 


NGC 7771 


4.0±0.3 


<2 
1.97 + 0.06 


1+1 


<2 


<1 


<2 



a Measured from high-resolution spectrum. 



Table 4. Spectral classification and other observables. 



Object 


8-13/tm 


Optical 






Infrared Photometry (Jy) 










-^IR 
M(H 2 ) 




AT 

N ame 


Class 


Class 


Ret 


12/tm 


25/tm 


60/tm 


100/tm 


Ret 


log 


L 


^ M / 


Ret 


Mkn 551 


PAH 


Hll 


3 


0.2415 


0.7986 


4.659 


5.686 


2 


11 


72 






IC 1632 


PAH 


Hll 


1 


0.98 


3.43 


22.19 


30.32 


1 


11 


54 


15.2 


1 


NGC 695 


PAH 


Hll 


1 


0.49 


0.8f 


7.61 


13.80 


1 


11 


64 


20.0 


1 


NGC 828 


PAH 


L 


6 


0.78 


f .04 


12.03 


24.38 


1 


11 


34 


15.1 


1 


MCG +05 — 06 — 036 


PAH 


Svl 

l_>J _L 


2 


0.2451 


.6944 


6.475 


11.45 


2 


11 


71 


10.1 


5 


IRAS 02528+4350 


Flat 






0.5046 


2.293 


2.762 




2 


12 


54 






IRAS 04154+1755 


PAH 


Sv2 


7 


0.2030 


0.7f4f 


3.822 


5.84 


2 


11 


75 






IRAS 05189 — 2524 


PAH 


Sv2 


1 


0.74 


3.50 


13.95 


12.52 


1 


12 


09 


54.3 


1 


NGC 2388 


PAH 


Hll 


1 


0.49 


f .96 


16.36 


29.68 


1 


11 


21 


21.4 


4 


NGC 2623 


PAH 


L or Sy2 


1 


0.24 


1.85 


25.72 


27.36 


f 


11 


55 


59.0 


1 


IRAS 09111-1007 


PAH 






0.17 a 


0.57 


7.19 


11.04 


f 


11 


99 


39 


f 


UGC 04881 


PAH? 


Hn 


1 


0.16 


0.62 


6.53 


10.21 


f 


11 


70 


17.6 


f 


IRAS 09433-1531 


PAH 






0.1876 


0.7174 


4.463 


4.698 


2 


11 


72 






NGC 3110 


PAH 


Hn 


1 


0.58 


1.10 


11.68 


23.16 


f 


11 


22 


9.3 


f 


IRAS 10565+2448 


PAH 


Hn 


1 


0.21 


1.21 


12.53 


16.06 


f 


11 


99 


48.3 


f 


Arp 299C 


PAH 






0.39 


2.2 






3 










Arp 299B1 


PAH+Sil 


Hn 


3 


1.7 


6.9 






3 










Arp 299C' 


PAH 






0.23 








4 










Arp 299 


PAH+Sil 






3.90 


24.14 


121.64 


122.45 


f 


11 


90 


71.7 


f 


Arp 299A 


PAH+Sil 


Hn 


3 


1.11 


12.4 


53 


64 


3,5 


11 


7 


125 


7 



Table 4 — Continued 



Object 


8-13/tm 


Optical 






Infrared Photometry (Jy) 








Mr 
M(H 2 ) 




Name 


Class 


Class 


Ref 


12/tm 


25/tm 


60/tm 


100/tm 


Ref 


w 


(M. ) 

V M / 


Ref 


UGC 08387 


PAH 


L 


1 


0.26 


1.37 


13.69 


24.90 


1 


11.58 


21 


3 


IRAS 13349+2438 


Flat 


Syl 


4 


0.6308 


0.8406 


0.6114 




2 


12.21 






UGC 08696 


Sil 


L 


1 


0.23 


2.30 


22.09 


22.44 


1 


12.14 


74.7 


1 


I Zw 107 


PAH 


Hn L? 


1 


0.22 


1.40 


9.15 


10.04 


1 


11.85 


40.6 


1 


NGC 6240 


PAH 


L 


1 


0.66 


3.67 


22.80 


27.35 


1 


11.83 


33.7 


1 


IRAS 17138-1017 


PAH 


Hn 


1 


0.62 


2.09 


15.84 


20.88 


1 


11.39 






IRAS 17208-0014 


PAH 


Hn 


1 


0.17 a 


1.71 


34.67 


37.65 


1 


12.40 


51 


2 


IRAS 18293-3413 


PAH 


Hn 


1 


1.01 


4.03 


35.92 


50.71 


1 


11.74 


15 


2 


ESO 593-IG008 


PAH 


L 


1 


0.20 


0.47 


6.59 


10.26 


1 


11.86 


22 


2 


IRAS 20210+1121 


Flat 


Sy2 


5 


0.2948 


1.398 


3.385 


2.682 


2 


11.87 


181 


6 


II Zw 96 


PAH? 


Hn 


1 


0.27 


2.28 


13.24 


11.61 


1 


11.89 






Zw 475.056 


PAH 


Sy2 


1 


0.32 


1.88 


8.75 


11.64 


1 


11.53 


37.3 


1 


NGC 7771 


PAH 


Hn 


1 


0.87 


2.18 


20.46 


37.42 


1 


11.35 


25.4 


1 



a From IRAS Faint Source Catalog (Moshir et al. 1990). 

References. — Optical spectral classifications are from (1) Veilleux et al. (1995), and (2) Mazzarella & Balzano (1986), (3) Mazzarella & 
Boroson (1993), (4) Beichman et al. (1986), (5) Perez et al. (1990), (6) Keel & Windhorst (1991), and (7) Crawford et al. (1996). Infrared fluxes 
are from (1) Soifer et al. (1989), (2) Moshir et al. (1990), (3) Wynn- Williams et al. (1991), (4) Keto et al. (1997), and (5) Joy et al. (1989). L m 
to M(H 2 ) ratios are from (1) Sanders, Scoville & Soifer (1991), (2) Mirabel et al. (1990), (3) Solomon et al. (1997), (4) Solomon & Sage (1988), 
(5) Kriigel, Chini & Steppe (1990), (6) Horellou et al. (1995), and (7) Sargent & Scoville (1991). 



Table 5. 60 /tin Flux limited subsample. 



Object 


8—13 /im 








Name 


Class 


Ref 


logim 


TTrr 1 sn^s fi\/TVn 9^1 ^ 

UL/Vj oUOo llvlKll zoi 1 


Sll 


A 


1 9 £Q 
1 Z.OO 


1 9£i 
— 1 .AO 


ir as 1 79ns nni a 


PAW 
r An 


1 


1 9 /in 


9 1 
— A . 1 


npri qqi / a ^t-, 9 9n^ 
ULv^j yyio 1 Arp zzui 


Sll 


A 


1010 
1 A . lo 


9 9 
— Z .Z 


TTPP SfiQfi fMVn 97^ 
U^Vj oOyO llvlKll ZlOl 


Sll 


1 


19 1/1 


9 n 
— z .u 


IP A C n^ICQ 9^9/1 


PA 11 

1 /ill 


1 


1 9 no 
1 z.uy 


1 Q 
— 1.0 




PA 11 
rAn 


1 


1 1 QQ 

1 1 .yy 


1 c 
— 1.0 


A m 9QQ 


pa 11 _l 

rAn -f- sil 


1 


I 1 Qfl 

I I .yu 


1 A Q 

— 1 .4y 


11 ZjW yo 


PA 11? 


1 


I 1 OQ 

I I .oy 


1 7 
— 1 . t 


INVj^ DZ4U 


PA 11 
rAn 


9 1 
Z , 1 


1 1 8^ 
1 1 .00 


1 ^ 

— 1.0 


IP A Q 1 C9QQ Q/11 Q 


PA 11 
r AH 


1 


1 1 7/1 

1 1 . 1 4 


1 

— 1.00 


IN 1 4oy 


P A 11 
r AH 


A 


1 1 fin 

1 1 .DU 


1 9/1 
— 1 .Z4 


T jpp 0007 
UUkj 000/ 


PA 11 
1 AH 


1 


1 1 

1 1 .00 


1 7 
— 1 . t 


W\j\^j 1D14 


P A 11 
1 AH 




1 1 

1 1 .00 


1 Qfi 
— 1 .OO 


NGC 2623 


PAH 


1 


11.55 


-2.0 


IC 1623 


PAH 


1 


11.54 


-1.4 


III Zw 35 


unknown 




11.54 


< -2.0 


NGC 34 


unknown 




11.41 


-1.6 


Mkn 331 


PAH 


2 


11.41 


-1.5 


ESO 507-G070 


unknown 




11.41 


-1.7 


IRAS 17138-1017 


PAH 


1 


11.39 


-1.4 


NGC 7771 


PAH 


1 


11.35 


-1.4 


NGC 828 


PAH 


1 


11.34 


-1.2 


NGC 1068 


sil 


2 


11.27 


-0.702 


NGC 3110 


PAH 


1 


11.22 


-1.3 


NGC 2388 


PAH 


1 


11.21 


-1.5 



References. — 



Classifications are from (l) this work or (2) RASW 



Table 6. Previously published results. 



Object 


EW(11.3) 


-^IR 
M(H 2 ) 




Name 


(/tm ) 


t L ® \ 

\ M @ I 


itel 


NGC 7469 


0.19 


46.6 


1 


NGC 1614 


0.36 


38.3 


1 


NGC 7714 


0.25 


23.0 


1 


NGC 0660 


0.52 


17. 


2 


NGC 6946 


0.49 


3.3 


2 


IC 342 


0.61 


3.0 


2 


NGC 0253 


0.24 


15 


2 


NGC 6240 


0.29 


33.7 


1 


IC 4553 (Arp 220) 


0.23 


81.9 


1 


Mkn 331 


0.70 


20.0 


1 


NGC 4194 


0.20 


66.1 


1 


NGC 2782 


0.54 


6.7 


1 



References. — Lib. to M(fi2) ratios are from (1) 
Sanders et al. (1991) or (2) Mirabel et al. (1990). 



